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INTRODUCTION 


In  a  previous  series  of  papers  [1-5],  matrix  methods  for  linear 
dynamic  analyses  of  lattice  structures  are  developed.  A  lattice  structure, 
in  this  context,  is  defined  to  be  an  idealized  network  of  one-dimensional 
members  which  are  connected  by  joints.  In  this  paper,  transfer  matrices 
and  joint  coupling  matrices  are  used  to  compute  the  natural  frequencies  of 
vibration  of  a  five-bay  planar  lattice  structure.  The  method  of  analysis 
is  applicable  to  general  two  and  three-dimensional  lattices.  The  necessary 
numerical  computations  may  be  performed  easily  using  a  personal  computer. 
Numerical  results  for  the  first  twenty-five  nonzero  natural  frequencies  of 
the  five-bay  lattice  structure  are  given  for  the  case  when  the  members  of 
the  lattice  are  modeled  as  Bernoulli-Euler  beams,  and  for  the  case  when 
the  lattice  members  are  modeled  as  Timoshenko  beams.  The  results  obtained 
here  are  compared  with  the  results  of  a  previous  analysis  [6]  using  a 
finite  element  method  and  an  experimental  modal  analysis.  A  short  dis¬ 
cussion  of  the  results  and  of  some  potential  applications  of  the  type  of 
analysis  presented  here  is  given. 


ANALYSIS 


Lattice  Structure  and  Model 

The  lattice  structure  considered  here  is  shown  in  Fig.  1.  The  lattice 
is  machined  from  a  single  piece  of  aluminum,  and  contains  no  welds  or 
fasteners. 

The  model  used  here  to  analyze  the  lattice  structure  of  Fig.  1  is 
shown  in  Fig.  2.  The  structure  is  modeled  as  an  idealized  lattice  of  one¬ 
dimensional  members  which  are  connected  by  joints.  It  is  assumed  that  the 
members  and  joints  can  move  only  in  the  plane  of  the  structure,  and  that 
all  motions  are  small. 

The  joints  of  the  lattice  model  are  labeled  1  through  12  as  shown  in 
Fig.  2.  It  is  assumed  that  each  joint  is  rigid,  and  that  each  connection 
between  a  joint  and  a  beam  is  rigid.  It  is  also  assumed  that  each  joint 
is  massless  and  has  no  spatial  extent.  The  assumptions  that  the  joints  are 
massless  and  have  no  spatial  extent  are  made  only  for  convenience;  some 
comments  about  the  analysis  of  joints  with  mass  and/or  spatial  extent  are 
given  below. 

It  is  assumed  that  all  members  of  the  lattice  model  of  Fig.  2  are 
identical,  and  that  each  member  can  extend  (and  contract)  axially  and  bend 
flexurally.  It  is  also  assumed  that  the  axial  and  flexural  motions  are 
uncoupled.  Two  different  member  models  are  used.  In  the  first  model, 
hereafter  called  the  Bernoulli-Euler  beam  model,  the  lattice  members  are 
modeled  as  classical  longitudinal  rods  for  axial  motions  and  as  Bernoulli- 
Euler  beams  for  flexural  motions.  In  the  second  model,  hereafter  called 
the  Timoshenko  beam  model,  the  lattice  members  are  modeled  as  classical 
longitudinal  rods  for  axial  motions  and  as  Timoshenko  beams  for  flexural 


motions.  In  both  models,  the  state  vector  at  any  point  x  of  a  lattice 


member  is  of  the  form 


z(x,t) 


(1) 


where  u(x,t)  is  the  longitudinal  displacement  of  the  member,  v(x,t)  is  the 
transverse  displacement  of  the  member,  ^(x,t)  is  the  rotation  of  the  member, 
M(x, t)  is  the  bending  moment  in  the  member,  V(x,t)  is  the  shear  force  in 
the  member,  F(x,t)  is  the  axial  force  in  the  member,  x  is  a  spatial  coordi¬ 
nate  which  extends  along  the  length  of  the  member  and  t  is  time.  The  com¬ 
ponents  of  the  state  vector  and  the  sign  convention  adopted  here  for  the 
components  of  the  state  vector  are  shown  in  Fig.  3.  Local  coordinate 
directions  x^  (i  *  1,2,  ...,  16)  are  assigned  to  the  lattice  members  as 
shown  in  Fig.  2. 


Joint  Coupling  Matrix  Relationships 

The  Fourier  transforms  z^  and  z^  of  the  state  vectors  and  shown 
in  Fig.  2  are  related  by  an  equation  of  the  form  [3] 


where  B^(co)  is  the  joint  coupling  matrix  of  joint  1.  Eqn.  (2)  is  the  joint 
coupling  matrix  relationship  for  joint  1.  Joint  coupling  matrix  relation¬ 
ships  for  joints  2  through  12  can  be  written  in  a  similar  manner  as 


B2(co) 


VW)i  >  =  ° 


Vw)l  ^9  }  =  ° 


B5(u>)  J  z12  >  =  0 


*15  >  -  * 


I 

to. 


The  locations  of  the  state  vectors  z_  (i  =  1,2,  . ..,  32)  are  shown  in 
Fig.  2.  The  joint  coupling  matrices  _B^  (i  =  1,2,  ...»  12)  may  be  derived 
from  the  general  formulas  given  in  [3],  or  directly  from  first  principles, 
as  is  done  in  Appendix  A.  The  joint  coupling  matrices  are  written  as  a 
function  of  radian  frequency  id  because  the  elements  of  the  joint  coupling 
matrices  depend,  in  general,  on  frequency.  (Note,  however,  that  for  the 
rigid,  massless  joints  considered  here,  the  elements  of  the  joint  coupling 
matrices  derived  in  Appendix  A  are  independent  of  frequency.)  The  deriva¬ 
tions  of  the  joint  coupling  matrices  in  Appendix  A  are  based  on  the  assump¬ 
tion  that  each  joint  is  completely  unconstrained.  Joint  coupling  matrices 
_Bj,  B^,  and  .B^  are  6  x  12  matrices,  and  joint  coupling  matrices  B^ 
through  B1q  are  9  x  18  matrices .  The  right  hand  side  of  each  of  eqns.  (2) 
through  (13)  is  zero  because  it  is  assumed  that  there  are  no  external  force 
or  moments  applied  to  the  joints. 

Eqns.  (2)  through  (13)  can  be  combined  into  a  single  equation  of  the 

form 


V“>2c  ‘  2 


(14) 


where  B  is  a  global  joint  coupling  matrix  given  by 


0 


0 


The  global  joint  coupling  matrix  B^(ui)  is  a  96  x  192  matrix,  and  the  global 

VJ 

state  vector  is  a  192  x  1  matrix.  Eqn.  (16)  contains  all  the  information 
about  the  dynamics  of  the  joints  in  the  lattice  model  of  Fig.  2  and  all  the 
connectivity  information  (that  is,  information  about  which  members  are 
connected  to  which  joints),  but  contains  no  information  about  the  dynamics 
of  the  members  of  the  lattice. 


Transfer  Matrix  Relationships 

The  Fourier  transforms  z^  and  t*ie  state  vectors  and  z^  in 

Fig.  2  are  related  by  an  equation  of  the  form 


-  9  - 


where  T(uj)  is  the  6  x  6  transfer  matrix  of  the  member  connecting  joints  1 
and  2.  Eqn.  (17)  is  the  transfer  matrix  relationship  for  the  member  con¬ 
necting  joints  1  and  2.  The  transfer  matrix  for  the  Bernoulli-Euler  beam 
model  and  the  transfer  matrix  for  the  Timoshenko  beam  model  are  given  in 
Appendix  B.  The  transfer  matrix  is  written  as  a  function  of  radian 
frequency  U)  because  the  elements  of  T^  depend,  in  general,  on  frequency. 
Transfer  matrix  relationships  for  each  of  the  remaining  members  in  the 
lattice  model  of  Fig.  2  can  be  written  as 


-23 

"  l(u)^18 

(27) 

^27 

=  l(u)z 22 

(28) 

-26 

=  T(u>)£25 

(29) 

—29 

=  1(w)I24 

(30) 

-32 

=  i(a))£2g 

(31) 

^31 

=  i^)I30 

(32) 

The  transfer  matrices  in  eqns.  (1?)  through  (32)  are  identical  because  the 
members  in  Fig.  2  are  identical,  and  because  of  the  choice  of  the  local 
coordinate  directions.  (As  discussed  in  [1],  reversing  the  sense  of  the 
local  coordinate  direction  changes  the  transfer  matrix.) 

Eqns.  (17)  through  (30)  can  be  combined  into  a  single  equation  of  the 

form 


h  '  k Mkll 


where  the  global  state  vector  is  given  by  eqn.  (16),  T^(u))  is  a  92  x  196 


global  transfer  matrix  given  by 


The  matrix  J_  in  eqn.  (34)  is  the  6 x 6  identity  matrix.  The  vector 
contains  half  the  state  vectors  in  Z^.  Eqn.  (35)  contains  all  the 
information  about  the  dynamics  of  the  members  in  the  lattice  model  of 
Fig.  2,  but  contains  no  information  about  the  joint  dynamics  and  no 
connectivity  information. 


|N| 
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Determination  of  Natural  Frequencies 


Substitution  of  eqn.  (33)  into  eqn.  (1A)  gives 


V“>V“)A;/2  ‘  -2 


where  the  product  Bc(co)T^(co)  is  a  96  x  96  matrix.  For  each  value  of  to, 
eqn.  (36)  is  a  system  of  homogeneous  linear  equations  for  the  components 
of  the  state  vector  The  values  of  <o  for  which  a  nontrivial  solution 

for  Z^/2  may  exist  must  satisfy  the  equation 


det(B^  (io)T^(u)) )  =  0 


The  values  of  co  which  satisfy  eqn.  (37)  are  the  natural  frequencies  of  the 
lattice  model  of  Fig.  2. 

Because  the  global  transfer  matrix  contains  trigonometric  and 
hyperbolic  functions  of  to,  eqn.  (37)  is  a  nonlinear  transcendental  equation 
in  to.  Since  the  members  of  the  lattice  model  of  Fig.  2  are  modeled  as  con¬ 
tinuous  beams,  there  is  an  infinite  number  of  values  of  co  which  satisfy 
eqn.  (37). 


^ v'  v*" ■  /*,*\*>‘  '  ** « ■ 
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NUMERICAL  RESULTS 


Method  of  Computation 

In  order  to  solve  eqn.  (37)  for  w,  the  determinant  of  B  („  )!,(..)  can 

~~G  ~~Q 

be  plotted  numerically  as  a  function  of  u),  and  the  values  of  w  for  which 
the  determinant  equals  zero  can  be  determined  graphically.  Alternatively, 
a  root-finding  algorithm  such  as  the  bisection  method  or  the  secant  method 
[7]  can  be  used. 

Both  the  graphical  method  and  the  root-finding  algorithms  require  a 
numerical  evaluation  of  det(B^(oj)’I-(co)) .  Since  most  of  the  elements  of 

Cj  (j 

B^(oj)  and  T  (to)  are  zero,  it  is  very  inefficient  to  encode  B„(lo)  and  T  (u) 
—\j  O  ■“ o  0 

into  a  computer  program  directly  from  eqns.  (15)  and  (34).  The  matrices 

B-(oo)  and  T^(o))  are  simple  enough  that  the  product  B-(u))T„  (oj)  may  be  com- 

— G  — G 

puted  in  6  x  6  block  form  by  hand.  The  result  of  such  a  hand  computation 


shows  that  IS  (w)T„(u))  is  a  banded  matrix.  Thus  the  determinant  of  B„(u.')T_(u:) 

Vj  ~ u  — O  — O 

may  be  computed  efficiently  by  using  a  computer  algorithm  specifically 

designed  for  banded  matrices. 


The  PASCAL  language  computer  program  used  here  to  evaluate 
det(B^(w)TG(u)) )  as  a  function  of  w  is  listed  in  Appendix  C.  The  heart  of 
the  program  is  the  procedure  bandet,  which  is  based  on  a  Gaussian  elimination 
algorithm  designed  to  solve  a  system  of  linear  equations  with  a  banded 
coefficient  matrix.  A  discussion  of  this  Gaussian  elimination  algorithm  is 
given  in  [8].  Only  the  portions  of  the  bandet  procedure  given  in  [8]  which 
are  necessary  to  evaluate  the  determinant  of  the  banded  matrix  are  used  here. 

Numerical  results  for  the  first  twenty-five  nonzero  natural  frequencies 
of  the  lattice  model  for  the  case  when  the  lattice  members  are  modeled  as 


Bernoulli-Euler  beams  and  for  the  case  when  the  lattice  members  are  modeled 


-  15  - 


as  Timoshenko  beams  are  given  in  Table  1.  The  results  in  Table  1  are 
obtained  by  first  plotting  det(IL,(oj)T^(io) )  as  a  function  of  to  in  steps  of 
10  rad/sec  to  evaluate  the  approximate  location  of  the  zero  crossings,  and 
then  evaluating  det  (B,(oj)T(oj)  )  as  a  function  of  to  in  steps  of  1  rad/sec  in 

Cj  Vj 

the  neighborhood  of  a  zero  crossing.  The  value  of  the  natural  frequency  is 
then  taken  as  the  average  of  the  two  values  between  which  det (B^(to)Tr,(uj) ) 
changes  sign.  The  material  and  geometric  properties  of  the  lattice  members 
which  are  used  in  the  computations  are  given  in  Appendix  B.  A  typical  plot 
of  det  (B,(oj)1\,(oj) )  as  a  function  of  oo  is  shown  in  Fig.  4.  The  square  data 
points  represent  the  values  of  det  (l^(oo)T(oo) )  which  were  actually  computed, 
and  the  curve  is  a  smooth  (second  order  polynomial)  fit  to  the  data  points. 
An  interesting  section  of  the  function  det(B  (oo)T^(oo))  is  shown  in 

(j  — 0 

Fig.  5.  It  appears  from  Fig.  5  that  det(B^(co)^(^) )  takes  the  value  zero 

Vj  — Vj 

somewhere  near  oj  =  3250  rad/sec.  The  very  small  slope  near  lo  -  3250  rad/sec 
suggests  that  the  zero  near  oj  =  3250  rad/sec  may  be  a  multiple  zero.  (A 
multiple  zero  of  multiplicity  k  is  a  value  of  oo  for  which  det  (B^(o))'T,(ol)  )  ■ 
and  the  first  k-1  derivatives  of  det(B„(oo)l_(oj))  with  respect  to  oj  are  also 

“"VJ  Vj 

equal  to  zero.)  However,  an  evaluation  of  det (B^(oo)l_(oo) )  near  oj  =  3250 

Vj  Vj 

rad/sec  with  a  smaller  step  size  for  oo  (see  Fig.  6)  shows  that  there  are  in 
fact  five  distinct  zero  crossings  between  oo  =  3242  rad/sec  and  oo  =  3260 
rad/sec.  (Note  that  the  vertical  scale  of  Fig.  6  is  much  different  from 
the  vertical  scale  of  Fig.  5.) 

The  natural  frequencies  computed  with  the  Bernoulli-Euler  beam  lattice 
model  are  plotted  as  a  function  of  flexible  mode  number  in  Fig.  7.  Fig.  7 
also  shows  the  natural  frequencies  obtained  in  a  previous  analysis  [6]  using 
a  finite  element  method  and  an  experimental  modal  analysis.  It  is  seen 


puted  here  and  those  given  in  [6].  The  maximum  difference  between  the 
results  obtained  here  and  those  given  in  [6]  is  six  percent. 


DISCUSSION 


The  number  of  natural  frequencies  in  Table  1  was  chosen  arbitrarily. 

The  method  presented  here  can  be  used  to  compute  as  many  natural  frequencies 
as  desired.  For  large  u,  the  value  of  det (B..(d)T  (uj) )  may  be  larger  than 
the  largest  number  which  a  particular  computer  can  store,  but  this  problem 
may  be  overcome  by  introducing  scale  factors  into  the  calculation  of  the 
determinant  [8],  As  shown  by  the  discussion  concerning  Figs.  5  and  6,  some 
care  is  required  in  finding  the  natural  frequencies,  since  they  may  be  very 
closely  spaced  in  some  frequency  ranges. 

The  results  in  Table  1  show  that  rotary  inertia  and  shear  deformation 
of  the  lattice  members,  which  the  Timoshenko  beam  model  includes  and  which 
the  Bernoulli-Euler  beam  model  does  not,  have  very  little  effect  on  the 
natural  frequencies  of  the  lattice  model  for  the  frequency  range  considered 
here.  Note,  however,  that  the  natural  frequencies  predicted  by  the  Timo¬ 
shenko  beam  model  are  less  than  or  equal  to  the  frequencies  predicted  by 
the  Bernoulli-Euler  beam  model,  and  that  the  difference  between  the  two 
models  increases  with  increasing  frequency.  It  can  be  shown  analytically 
that  for  a  single  beam  which  is  simply  supported  at  each  end,  the  effect 
of  rotary  inertia  and  shear  deformation  is  to  decrease  the  natural  frequencies 
of  flexural  vibration,  and  that  the  effect  increases  with  increasing 
frequency  [9]. 

The  joint  models  adopted  here  are  perhaps  the  simplest  possible.  How¬ 
ever,  the  excellent  agreement  between  the  results  computed  here  and  the 
results  of  the  experimental  analysis  in  [6]  shows  that  the  simple  joint 
models  are  certainly  useful  in  the  frequency  range  considered  here.  Also, 
more  complicated  joint  models  may  be  included  in  the  present  analysis  with 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  type  of  analysis  presented  here  may  be  applied  to  any  system  which 
may  be  modeled  as  a  connected  network  of  one-dimensional  members .  (Note, 
for  example,  that  this  kind  of  analysis  is  not  restricted  to  structural  or 
mechanical  systems,  and  is  related  to  techniques  used  in  the  analysis  of 
electrical  and  microwave  networks  [10}.)  The  necessary  numerical  computa¬ 
tions  are  straightforward,  and  are  easily  performed  with  a  personal  computer. 

As  discussed  above,  the  effect  of  various  member  models  and  joint  models 
on  the  natural  frequencies  of  the  lattice  may  be  determined  by  simply 
altering  the  appropriate  transfer  matrices  and  joint  coupling  matrices. 

Also,  and  perhaps  more  importantly,  the  effects  of  changes  in  the  structure 
of  the  lattice  may  be  considered  using  the  techniques  presented  here. 

Changes  in  structure  may  be  due,  for  example,  to  disconnected  joints  or 
damaged  members.  Thus  the  techniques  discussed  here  may  be  used  to  begin 
to  study  the  following  nondestructive  evaluation  question:  given  a  mea¬ 
surement  of,  say,  certain  natural  frequencies  of  a  lattice  structure,  what 
is  it  possible  to  conclude  about  the  structural  integrity  of  the  lattice? 

This  nondestructive  evaluation  question  may  also  be  expressed  in  the  wider 
context  of  the  system  identification  problem,  in  which  the  properties  of 
an  unknown  system  are  deduced  from  certain  measurements  of  the  system. 
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TABLE  1  First  twenty-five  nonzero  natural  frequencies  of  lattice  model. 


Flexible 

Mode 

Number 


Bernoulli-Euler  Beams 
(rad/sec) 


Timoshenko  Beams 
(rad/sec ) 


308.5 

308.5 

433.5 

433.5 

606.5 

605.5 

726.5 

725.5 

932.5 

931.5 

1438.5 

1437.5 

1535.5 

1534.5 

1766.5 

1764.5 

2071.5 

2067.5 

2115.5 

2111.5 

2307.5 

2303.5 

2479.5 

2474.5 

2887.5 

2881.5 

2979.5 

2973.5 

3242.5 

3232.5 

3249.5 

3238.5 

3254.5 

3243.5 

3257.5 

3246.5 

3258.5 

3248.5 

3552.5 

3544.5 

4120.5 

4108.5 

5923.5 

5903.5 

6094.5 

6072.5 

6704.5 

1 

6675.5 

1  6894.5 

6861.5 

T  T 


Typical  plot  of  det  (B^ (cu)  T ^  (uj)  )  as  a  function  of  oo,  obtained  using  Bernoull 
Euler  beam  mode] . 


7.5  x  10 


obtained  using  Bernoull i-Eul e’r  beam  model. 


Frequency  (Hz) 


2 


•I* 


T 


T 


Fig.  7  Natural  frequencies  obtained  by  finite  element  analysis, 
experimental  modal  analysis,  and  Bernoull i-Euler  beam 
lattice  model. 
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APPENDIX  A:  DERIVATION  OF  JOINT  COUPLING  MATRICES 


In  this  appendix,  the  joint  coupling  matrices  for  the  lattice  model 
of  Fig.  A1  are  derived.  The  joints  of  the  lattice  model  are  labeled  1 
through  12  as  shown  in  Fig.  Al.  It  is  assumed  that  all  joints  and  members 
can  move  only  in  the  plane  of  the  lattice,  and  that  all  motions  are  small. 
It  is  also  assumed  that  each  joint  is  rigid,  massless,  and  has  no  spatial 
extent,  and  that  each  connection  between  a  member  and  a  joint  is  rigid. 

The  state  vector  at  any  point  x  of  any  member  of  the  lattice  model  of 
Fig.  Al  is  assumed  to  be  of  the  form 


z.(x,t) 


ru(x,tr 

i  v(x,t) 
|i|»(x,t) 
M(x,t) 
V(x, t) 
F(x,t) 


(Al) 


where  u(x,t)  is  the  longitudinal  displacement  of  the  member,  v(x,t)  is  the 
transverse  displacement  of  the  member,  ijj(x,t)  is  the  rotation  of  the  member 
M(x,t)  is  the  bending  moment  in  the  member,  V(x,t)  is  the  shear  force  in 
the  member,  F(x,t)  is  the  axial  force  in  the  member,  x  is  a  spatial  coor¬ 
dinate  which  extends  along  the  length  of  the  member  and  t  is  time.  The 
components  of  the  state  vector  and  the  sign  convention  adopted  here  for 
the  components  of  the  state  vector  are  shown  in  Fig.  A2.  Local  coordinate 
directions  x^  (i  =  1,2,  ...,  16)  are  assigned  to  the  lattice  members  as 
shown  in  Fig.  Al.  Throughout  this  appendix,  an  overbar  will  denote  a 


Fourier  transform. 


Joint  1 


The  components  of  the  state  vectors  and  z?  and  a  free-body  diagram 
of  joint  1  are  shown  in  Fig.  A3.  The  equilibrium  equations  for  joint  1 


M  +  M2  =  0 


V1  -  F2  "  ° 


F1  +  V2  -  0 


Eqns.  (A2),  (A3)  and  (A4)  can  be  written  in  matrix  form  as 


1  0  0  M 


0  0  1 


1  0  0  M 


0  1  0  I  v2  V  +  0  0-l  <V2V*<0 


J  L 


0  10 


The  compatibility  equations  for  joint  1  are 


u^  +  v2  =  0 


V1  -  u2  "  0 


*1  -  4'2  ■  0 


Eqns.  (A6),  (A7)  and  (A8)  can  be  written  in  matrix  form  as 


10  0 


0  1  0 


0  0  1 


0  10 


V1  f  +  i'1  0  ° 


0  0-1 


’2-1° 


Combining  eqns.  (A5)  and  (A0)  and  taking  the  Fourier  transform  of  the 


resulting  equation  give 


.'Tiv.'vW'  hjt 


10000001000  o' 
010000  -1  00000 
00100000  -1  000 
000100000100 
00001000000  -1 
000001000010 


(A10) 


The  6  x  12  matrix  in  eqn.  (A10)  is  the  joint  coupling  matrix  of  joint  1. 

Joint  2 

The  components  of  the  state  vectors  and  z_^  and  a  free-body  diagram 
of  joint  2  are  shown  in  Fig.  A4 .  The  equilibrium  equations  for  joint  2  are 


1  0  0  M, 


-10  0 


oio  jv3f+  001 


0  0  1 


0-10 


(All) 


The  compatibility  equations  for  joint  2  are 


•  ’JiV  W  V  "Vji  *  A  ’ji  '  * 


.  'Sv/y'-AA-A. 


.  V-. 'vVa'Iv'A 


M,  -  M,  -  M  =  0 
DO/ 

V5  -  v6  -  F?  =  0 

F5  -  F6  +  V7  =  ° 


( A 1 4  ) 

( A 1 5 ) 

(Alh) 


Eqns.  (A14) ,  (A15)  and  (A16)  can  be  written  in  matrix  form  as 


— 

1 

0 

0 

M5 

-1 

0 

~1 

0 

0 

1 

0 

V5 

■  + 

0 

-1 

0 

0 

0 

1 

.F5 

0 

0 

-1 

Mb' 

— 

-1 

0 

0 

V 

" 

0 

V6 

•  + 

0 

0 

-1 

V. 

/ 

•  =  ’ 

0 

.V 

0 

1 

0 

,F7. 

0 

(A17) 


The  compatibility  equations  for  joint  3  are 


U5  "  U6 


=  0 


V  —  v 

5  6 


=  0 


*5  -  *6 


=  0 


u5  ~  v7 


=  0 


v5  +  u7 


=  0 


**5  -  *7 


=  0 


(A18) 

(A19 ) 

(A20) 

(A21) 

(A22 ) 

( A  2  3 ) 


Eqns.  (A18)  through  (A23)  can  be  written  in  matrix  form  as 


resulting  equation  give 


I  1 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1  0 

1 

l 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

a 

0 

0 

0 

0 

0 

0 

t 

i  u 

0 

1 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

!  1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

!  u 

l 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1  0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

'  c 

i 

0 

0 

1 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

-1 

0 

0 

1  0 
\ 

0 

0 

0 

1 

0 

0 

0 

0 

0 

-1 

Q 

0 

0 

0 

0 

0 

-1 

i  0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

1 

0 

The  9  x  18  matrix  in  eqn.  (A2b)  is  the  joint  coupling  matrix  of  joint  3. 


Joint  4 

The  components  of  the  state  vectors  z_g ,  and  and  a  free-body 

diagram  of  joint  4  are  shown  in  Fig.  A6.  The  equilibrium  equations  for 
joint  4  are 


~  1 
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0~ 

M 

8 

~  i 
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O- 

V 

-1 

0 

0" 

M10 

o' 

0 

1 

0 

< 

V8 

+ 
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-1 

V9 
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V10 
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00 

0 

1 

0 

F9 
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—  X 

0 

F10 

0 

1.  < 

(A27) 


The  compatibility  equations  for  joint  4  are 


~  1 

0 

0 

00 

3 

"0 

1 

0 

u9 

0 

0 

1 

0 

00 

> 

.  + 

-1 

0 

0 

vn 
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.  =  i 

0 

0 

0 

1 

^8 

0 

0 

-1 

'  9 

0 

— 

— 

— 

— 

1 

0 

o- 

c 

00 

"0 

1 

0“ 

u10 

' 

0 

0 

1 

0 

00 

> 

,  + 

-1 

0 

0 

o 

t-H 

> 

*  =  i 

0 

_° 

0 

1 

,T8. 

0 

0 

-1 

.'10. 

0 

( A  2  8 ) 


( A  2  9 ) 


Combining  eqns.  (.A27),  (A28)  and  (A29)  and  taking  the  Fourier  transform 


of  the  resulting  equation  give 


# 


* 


The  9  x  18  matrix  in  eqn.  (A30)  is  the  joint  coupling  matrix  B,  of  joint  4. 


Joints  5,  7  and  9 

The  geometry  of  each  of  joints  5,  7  and  9  is  identical  to  the  geometry 
of  joint  3.  Also,  the  relative  directions  of  the  local  coordinate  direc¬ 
tions  of  the  members  attached  to  each  of  joints  5,  7  and  9  are  identical  to 
the  relative  directions  of  the  members  attached  to  joint  3.  Therefore, 
each  of  the  joint  coupling  matrices  EL ,  E5^  and  is  equal  to  B^. 

Joints  6,  8  and  10 

The  geometry  of  each  of  joints  6,  8  and  10  is  identical  to  the  geometry 
of  joint  4.  Also,  the  relative  directions  of  the  local  coordinate  direc¬ 
tions  of  the  members  attached  to  each  of  joints  6,  8  and  10  are  identical 
to  the  relative  directions  of  the  members  attached  to  joint  4.  Therefore, 
each  of  the  joint  coupling  matrices  B^,  Bg  and  _B^q  is  equal  to  B^. 


Joint  11 

The  components  of  the  state  vectors  z^,g  and  z^q  and  a  free-body  diagram 
of  joint  11  are  shown  in  Fig.  A7.  The  equilibrium  equations  for  joint  11 
are 


MBWWflwr  fygwwwc  ctwwmh  WWW!  wcototw 


The  compatibility  equations  for  joint  11  are 


1 

0 

0 

f  ' 

U29 

0 

-1 

0 

u30 

0 

0 

1 

0 

V29 

*  + 

1 

0 

0 

< 

v30 

,  =  < 

0 

0 

0 

1 

_ 

v29_ 

0 

0 

-1 

_  '^30, 

0 

.  * 

(A32) 


Combining  eqns.  (A31)  and  (A32)  and  taking  the  Fourier  transform  of  the 


resulting  equation  give 


10  0 
0  10 
0  0  1 
0  0  0 
0  0  0 
0  0  0 


0  0  0 
0  0  0 
0  0  0 
10  0 
0  10 
0  0  1 


0-10 
10  0 
0  0-1 

0  0  0 

0  0  0 

0  0  0 


0  0  0 

0  0  0 

0  0  0 

-10  0 

0  0-1 

0  10 


(A33) 


The  6  x  12  matrix  in  eqn.  (A33)  is  the  joint  coupling  matrix  B, ,  of  joint 


Joint  12 


The  components  of  the  state  vectors  z  ^  and  2^2  anci  a  free-body  diagram 
of  joint  12  are  shown  in  Fig.  A8.  The  equilibrium  equations  for  joint  12 
are 
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The  compatibility  equations  for  joint  12  are 


1 

0 

0 

U31 

0 
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,  + 

-1 
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0 

V32 
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-1. 
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(A3  A ) 


(A35) 


Combining  eqns.  (A3A)  and  (A35)  and  taking  the  Fourier  transform  of 


the  resulting  equation  give 


1000000100 
010000  -1  000 
00100000  -1  0 
0001000001 
0000100000 
0000010000 


0  0 

0  0 

0  0 

0  0 

0  -1 

1  0 


(A36 


6  x  12  matrix  in  eqn.  (A36)  is  the  joint  coupling  matrix  _B  9  of  joint 


model  . 
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APPENDIX  B:  TRANSFER  MATRIX  RELATIONSHIPS 

In  this  appendix,  the  transfer  matrix  relationships  for  two  lattice 
member  models  are  given.  In  the  first  model,  hereafter  called  the  Bernoulli 
Euler  beam  model,  the  lattice  member  is  modeled  as  a  classical  longitudinal 
rod  for  axial  motions  and  as  a  Bernoull i-Euler  beam  for  flexural  motions. 

In  the  second  model,  hereafter  called  the  Timoshenko  beam  model,  the  lattice 
member  is  modeled  as  a  classical  longitudinal  rod  for  axial  motions  and  as 
a  Timoshenko  beam  for  flexural  motions.  It  is  assumed  in  both  models  that 
the  axial  and  flexural  motions  are  uncoupled.  For  both  models,  the  state 
vector  at  any  point  x  of  the  lattice  member  is  of  the  form 


z(x,t)  = 


(Bl) 


where  u(x,t)  is  the  longitudinal  displacement  of  the  member,  v(x,t)  is  the 
transverse  displacement  of  the  member,  +(x,t)  is  the  rotation  of  the  member, 
M(x,t)  is  the  bending  moment  in  the  member,  V(x,t)  is  the  shear  force  in 
the  member,  F(x,t)  is  the  axial  force  in  the  member,  x  is  a  spatial  coordi¬ 
nate  which  extends  along  the  length  of  the  member  and  t  is  time.  The 
components  of  the  state  vector  and  the  sign  convention  adopted  here  for 
the  components  of  the  state  vector  are  shown  in  Fig.  Bl .  In  Fig.  Bl ,  the 
left  end  of  the  member  is  designated  as  point  1,  and  the  right  end  of  the 


member  is  designated  as  point  2.  The  local  coordinate  x  is  defined  as 
shown  in  Fig.  Bl.  Throughout  this  appendix,  an  overbar  will  denote  a 
Fourier  transform.  Derivations  of  the  transfer  matrix  relationships  given 
here  can  be  found  in  [1]  and  [11]. 

Bernoulli-Euler  Beam  Model 

The  transfer  matrix  relationship  for  the  Bernoulli-Euler  beam  model 
is  given  by 


where 


In  eqns .  (B2)  through  (B8) ,  u  is  the  mass  per  unit  length  of  the  member, 
u>  is  radian  frequency,  Z  is  the  length  of  the  member  from  point  1  to  point 
2,  E  is  the  elastic  modulus  of  the  member,  El  is  the  flexural  rigidity  of 
the  member,  A  is  the  cross-sectional  area  of  the  member  and  p  is  the  mass 
density  (mass  per  unit  volume)  of  the  member.  For  numerical  computations, 
the  following  values  of  the  material  and  geometric  constants  are  used: 

U  =  2.44  x  lCf5  lb-sec2/in2,  l  =  9.85  in,  E  =  10  x  106  lb/in2,  El  =  4880 
lb-in  ,  A  =  9.38  x  10  2  in2  and  p=  2.6  x  10  ^  lb-sec2/in^. 


Timoshenko  Beam  Model 

The  transfer  matrix  relationship  for  the  Timoshenko  beam  model  is 
given  by 


54 


In  eqns.  (B9)  through  (B20) ,  y  is  the  mass  per  unit  length  of  the  member, 
w  is  radian  frequency,  £  is  the  length  of  the  member  from  point  1  to  point 


2 ,  G  is  the  shear  modulus  of  the  member,  A  is  the  cross-sectional  area 

s 

of  the  member  divided  by  the  geometric  correction  factor  K,  i  is  the 
radius  of  gyration  of  the  member  cross  section,  El  is  the  flexural 
rigidity  of  the  member,  p  is  the  mass  density  (mass  per  unit  volume)  of 
the  member,  E  is  the  elastic  modulus  of  the  member  and  A  is  the  cross- 
sectional  area  of  the  member.  For  numerical  computations,  the  following 
values  of  the  material  and  geometric  constants  are  used:  y  =  2.44  x  10  ^ 
lb-sec2/in2,  £  =  9.85  in,  G  =  3.8  x  106  lb/in2,  Ag  =  0.113  in2,  <  =  0.833, 
i  =  7.21  x  10  2  in,  El  =  4880  lb-in2,  p  =  2.6  x  10  ^  lb-sec2/in^, 

E  =  10  x  10^  lb/in2  and  A  =  9.38  x  10  2  in2. 


PROGRAM  detcalc  ( input , output) ; 


(*  This  program  is  used  to  evaluate  de t ( BG ( w ) TG ( w )  )  for  a  *) 
(*  given  value  of  w.  The  nonzero  elements  of  BG(w)TG(w)  are  *) 
(*  stored  in  the  matrix  STBGTG.  The  matrix  DIAG  contains  *) 
(*  the  nonzero  diagonal  lines  of  BG(w)TG(w).  The  matrix  *) 
(*  DIAG  is  passed  to  the  procedure  bandet,  which  returns  *) 
(*  the  value  of  de t  (  BG (  w)  TG  ( w) )  .  A  discussion  of  the  *) 
(*  procedure  bandet  is  given  in  Linear  Algebra,  by  J.  H.  *) 
(*  Wilkinson  and  C.  Reinsch,  Springer-Verlag ,  1971.  * ) 


CONST 

m  =  6  ; 
n  =  96  ; 
p  =  36 ; 
ml  =23; 
m2  =17; 

E  =  10E6; 

G  =  3-8E6; 

1th  =  9.85; 
ro  =  2.6E-4; 
kap  =  0.8333; 

El  =  4880; 
ar  =  0.0938; 

AS  =  0.1126; 
mu  =  2.439E-5; 
r  =  0.0721 ; 

TYPE 

mmatrix  =  ARRAY  [l..m,1..m]  OF  real; 
npmatrix  =  ARRAY  [l..n,1..p]  OF  real; 
dmatrix  =  ARRAY  [ 1 . . n  ,-ml . . m2  ]  OF  real; 

VAR 

th,w,b4,sig,tau,lam1 ,lam2,lam,c0,c1 ,c2,c3,det,b 
STBGTG 
DIAG 

T,ID,B12,B22,B31 , B32 , B33 , B34 , B35 , B36 , B41 ,B42, 
B45,B46,B112,B122 
i  »  j  i  k 

Timoshenko,  BernoulliEuler 
PROCEDURE  bandet( n ,ra1 , m2 : integer ;  VAR  A:dmatrix;  VAR  detrreal); 
VAR 

dl , x , norm , macheps  :real; 
i  ,  j , k , 1  : integer ; 


:  rea  1 ; 
.•npmatrix ; 
:  dmatrix ; 

smmatrix ; 

:  integer ; 

: Boolean ; 
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BEGIN 

macheps 

norm 


:=  0.000000000000001; 

:=  0.0; 


FOR  i  :=  1  to  n  DO 
BEGIN  x  :=  0.0; 

FOR  j  :=  -ml  to  m2  DO 
x  :=  x  +  abs ( A [ i , j ] )  ; 

IF  (norm  <  x)  THEN  norm  :=  x; 

END; 

I  :=  ml  ; 

FOR  i  :=  1  to  ml  DO 
BEGIN 

FOR  j  :=  1-i  to  m2  DO 
A[i,  j-1]  :=  A[i, j] ; 

1  :=  1  -  1  ; 

FOR  j  :=  m2-l  to  m2  DO 
A [ i »  j  ]  8  =  0.0; 

END; 

dl  s=  1.0;  1  :=  ml; 

FOR  k  :=  1  to  n  DO 

BEGIN  x  :=  A[k,-ml];  i  :=  k; 

IF  ( l<n )  THEN  1  :=  1  +  1; 

FOR  j  :=  k  +  1  to  1  do 
IF  abs ( A[ j ,-ml ] )  >  abs(x)  THEN 
BEGIN  x  :=  A [ j , -ml ] ;  i  : *  j;  END; 
dl  :=  dl  *  x; 

IF  (x  =  0)  THEN 
BEGIN 

A[k,-ml]  :=  norm  *  macheps; 

END  ; 

IF  (i  <>  k)  THEN 
BEGIN  dl  :=  -dl ; 

FOR  j  :=  -ml  to  m2  DO 
BEGIN  x  :=  A[k,j];  A[k,j]  :=  A[i,j] 
A [  i  ,  j  ]  j=  x  ; 

END; 

END; 

FOR  i  :=  k  +  1  to  1  DO 
BEGIN  x  :=  A [ i , -ml ] / A [ k , -ml  ]  ; 

FOR  j  :=  1-ml  to  m2  DO 
A [ i *  0-1 J  :=  A[ i  ,  j  ]  -  (x  *  A [ k , j  ]  )  ; 

A [  i  , m2  ]  8=  0  ; 

END; 

END; 

det  s  =  dl  ; 

END; 

FUNCTION  sinh  (xsreal)  :  real; 

BEGIN 

sinh  :=  ( exp ( x )  -  exp(-x))/2.0 
END; 
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FUNCTION  cosh  (x:real)  :  real; 
BEGIN 

cosh  :  =  (exp(x)  +  exp(-x))/2.0 
END; 

BEGIN 

Timoshenko  :=  false; 

BernoulliEuler  :  =  true; 


FOR  i  :=  1  to  m  DO 
BEGIN 

FOR  j  :=  1  to  m 
BEGIN 

ID[i,j]  : 
B12[i,j]  : 

B22[i,j]  : 
B3l[i,j]  : 
B36[i,j]  : 
B42[i,j]  : 
B46[i,j]  : 
B1 22 [ i  ,  j  ]  : 
END; 


(*  initialize  submatrices  of  BGTG  *) 


B32[i,  j] 
B33 [ i » 0  ] 
B34[i.  J] 
B35 [ i » 0  ] 
B41 [i,  J] 
B45[i,  j] 
B112[i,  j] 


=  0.0; 

=  0.0; 

=  0.0; 

=  0.0; 

=  0.0; 

=  0.0; 

:  =  0.0; 


(*  form  the  joint  coupling  matrices  *) 


B12[1 ,2] 

=  1  .0 

B1  2  [  4  1  4  ] 

=  1  .0 

B1  2  [  2 , 1 

=  -1  .0 

B1 2  [  5  »  6  ] 

=  -1.0 

B12 [3,3] 

=  -1  .0 

B1 2 [6 , 5  ] 

=  1  .0 

B22  1 , 2 J 

=  1  .0 

B22 [ 4  >  4  J 

=  -1  .0 

B22 [2,1 

=  -1  .0 

B22  [5,6] 

=  1.0 

B22[3,3] 

=  -1  .0 

B22 [6,5] 

=  -1  .0 

B31 [1  ,1  J 

=  1  .0 

B31 [4,1 

=  1 .0 

B31 [2,2] 

=  1  .0 

B31 [5,2] 

=  1  .0 

B31 [3,3] 

=  1  .0 

B31 [6,3] 

=  1  .0 

B32[1,1] 

=  -1  .0 

B33 [ 4  >  2  ] 

=  -1.0 

B32 [ 2 , 2  ] 

=  -1  .0 

B33 [5  » 1  ] 

=  1  .0 

B32 [ 3 , 3 ] 

=  -1  .0 

B33 [6 , 3  ] 

=  -1  .0 

B34[1 ,4] 

=  1  .0 

B35 [1,4] 

=  -1  .0 

B34 [2,5] 

=  1  .0 

B35 [2,5] 

=  -1.0 

B34 [3,6] 

=  1.0 

B35 [3,6] 

=  -1.0 

B36[1 ,4] 

=  -1  .0 

B41 [4,1 

=  1  .0 

B36[ 2,6] 

=  -1.0 

B4I [5,2] 

=  1  .0 

B36 [3,5] 

=  1  .0 

B41  [6,3] 

=  1.0 

B42 [ 4 1 2  ] 

=  1  .0 

B45 [4,4] 

=  1  .0 

B42 [5,1] 

=  -1  .0 

B45 [5,6] 

=  -1  .0 

B42 [6,3] 

=  -1.0 

B45 [ 6 , 5  ] 

=  1  .0 

B46  1 ,2] 

=  1  .0 

B46 [4,4] 

=  -1  .0 

B46 [2,1] 

=  -1  .0 

B46 [5,6] 

=  1  .0 

B46 [3,3] 

=  -1  .0 

B46 [6,5] 

=  -1  .0 
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B112[1 ,2] 

=  -1.0;  B112[4.4J 

=  -1.0 

B1 1 2  [  2 , 1  ] 

=  1.0;  B1 1 2 [ 5  »  6 ] 

=  -1  .0 

B1 1 2  [  3  *  3  J 

=  -1.0;  B1 1 2 [6 , 5  J 

=  1  .0 

B122[1  ,2] 

=  1.0;  B122[4.4] 

=  1  .0 

B1 22 [ 2 , 1 J 

=  -1.0;  B122[5,6] 

=  -1.0 

B1 22 [ 3  >  3  ] 

=  -1.0;  B1 22 [6 , 5 ] 

=  1  .0 

ID[1,1  J 

=  1.0;  ID [ 4 , 4 ] 

=  1.0 

IDL  2  » 2 ] 

=  1.0;  ID [ 5  *  5 ] 

=  1 .0 

ID[ 3 » 3  ] 

=  1.0;  ID[6 , 6] 

=  1  .0 

w  :=  10.0;  (*  assign  a  value  to  w  *) 

WHILE  (w  <  10000.0)  DO 
BEGIN 

FOR  i  :=  1  to  n  DO  (*  initialize  the  matrix  STBGTG  *) 

BEGIN 

FOR  j  :=  1  to  p  DO 
BEGIN 

STBGTG [  i  ,  j  ]  :=  0.0; 

END; 

END; 

IF  (Timoshenko)  THEN 

BEGIN  (*  form  transfer  matrix  of  Timoshenko  beam  *) 
th  :=  1th  *  w  *  sqrt(ro/E); 
b4  j  =  mu  *  sqr(w)  *  sqr ( sqr ( 1th) ) /El ; 

sig  :=  mu  *  sqr(w)*  sqr ( 1th) / ( G*as ) ; 
tau  :=  mu  *  sqr(r)*  sqr(w)*  sqr(lth)/EI; 
laml  :=  Bqrt(sqrt(b4  +  ( ( sqr ( Big- tau ) ) /4 ) ) 

-  (sig+tau)/2) ; 

lam2  :=  sqrt(sqrt(b4  +  ( ( sqr ( sig- tau) ) /4 ) ) 

+  ( Big+tau)/ 2) ; 

lam  :=  l/(sqr(laml)  +  sqr(lam2)); 
cO  :=  lam  *  (sqr(lam2)  *  cosh(  laml ) 

+  sqr(laml)  *  cos(lam2)); 
cl  :=  lam  *  ( sqr ( lam2 )/ laml  *  sinh(laml) 

+  sqr ( laml )/ lam2*  sin(lam2)); 
c2  :=  lam  *  (cosh(laml)  -  cos(lam2)); 
c3  :  =  lam  *  (1/laml  *  sinh(laml) 

-  l/lam2  *  sin(lam2)); 

T[1 ,1  ]  :=  co  b ( th ) ; 

T [ 1  ,2]  :=  0.0; 

T[1 ,3]  :=  0.0; 

T [ 1  ,4]  s=  0.0; 

T[1,5]  :=  0.0; 

T [ 1 »  6 ]  :=  1th/ ( E*ar ) *sin ( th ) / th ; 

T [ 2 , 1 J  :=  0.0; 

T [ 2 , 2  ]  :=  cO  -  ( sig*c2 ) ; 

T [ 2 , 3  ]  :=  lth*(c1  -  ((sig  +  tau)  *  c3)); 

T  [  2 , 4  ]  :=  (sqr(lth)/EI)*c2; 


T[2,5J  :=  (sqr(lth)*lth/( b4*EI ))*(-( sig* cl  ) 

+  ( ( b4  +  sqr(sig) )*c3) ) ; 

T[ 2 , 6]  :=  0.0; 

T  3 , 1  s  =  0.0; 

T[ 3 » 2 ]  :=  (b4/lth)*c3; 

Tf 3 > 3 ]  :=  c0  -  ( tau*c2 ) ; 

T .3,4.  :=  (lth/El)  *  (cl  -  (tau*c3)); 

T[ 3 , 5 ]  :=  T [ 2 , 4 ]  ; 

T [ 3  » 6 ]  :=  0.0; 

T[ 4  » 1 ]  :=  0.0; 

T[ 4  » 2 ]  :=  b4  *  ( El/sqr ( 1th) )  *  c2; 

T  4»3  :  =  (El/lth)  *  (-(tau*cl)  +  ( ( b4  +  sqr ( tau ) ) * c3 ) ) ; 

T[4.4]  :=  T[3,3j; 

T[ 4 i 5 ]  :=  T[2,3]  ; 

T[4,6]  :=  0.0; 

T  5 , 1  : =  0.0; 

T [ 5  *  2 ]  j=  b4*(El/(sqr(lth)*lth))*(c1  -  (sig*c3)); 

T[5,3]  :=  T[ 4  » 2 ] ; 

T[ 5  » 4  ]  :=  T [ 3  » 2  ]  ; 

T[ 5  »  5  ]  :=  T [ 2 , 2  J  ; 

T[5,6J  :=  0.0; 

T[6,l]  j=  -xnu*lth*sqr  (  w)  *sin  ( th)/ th  ; 

T [ 6 , 2 ]  :=  0.0; 

T[6 , 3 ]  :=  0.0; 

T[6,4]  :=  0.0; 

T [6 , 5 ]  :=  0.0; 

T[6 ,6  :=  cos( th) ; 

END; 

IF  (BernoulliEuler)  THEN 

BEGIN  (*  form  transfer  matrix  of  Bernoulli-Euler  beam  *) 
b4  :=  mu  *  sqr(w)  *  sqr ( sqr ( 1th) )  /  El; 

b  :=  sqrt ( sqrt ( b4 ) ) ; 

th  :=  1th  *  w  *  sqrt(ro/l); 

cO  :=  0.5  *  (cosh(b)  +  cos(b)) ; 

cl  :=  (1  /  (2  *  b))  *  (sinh(b)  +  sin(b)); 

c2  :=  (1  /  (2  *  sqr(b)))  *  (coBh(b)  -  cos(b)); 

c3  s*  (1  /  (2  *  sqr(b)  *  b))  *  (sinh(b)  -  sin(b)); 

T [ 1 , 1  ]  :=  cos ( th) ; 

T[1 ,2]  :=  0.0; 

T[1,3]  :=  0.0; 

T [ 1 » 4 ]  :=  0.0; 

T[1,5]  :=  0.0; 

T [ 1  * 6 ]  :=  1th  /  (E  *  ar)  *  sin(th)  /  th ; 

T [ 2 , 1  ]  :=  0.0; 

T [ 2 , 2 ]  :=  cO; 

T[ 2 , 3  ]  :=  1th  *  cl  ; 

T[ 2 , 4 ]  :=  sqr ( 1th )  *  c2  /  El; 

T[ 2 , 5 ]  :=  sqr(lth)  *  1th  *  c3  /  El; 

T[2,6]  :=  0.0; 

T[ 3 , 1 ]  :=  0.0; 


T[ 3 » 2 ]  :=  b4  *  C3  /  1th; 

T[3,3]  :=  CO; 

T [ 3  *  4  3  :=  1th  *  cl  /  El; 

T[ 3  »  5 ]  :=  T[ 2 ,4  ]  5 
T[ 3*6]  :=  O.O; 

T[ 4  *  1 ]  :=  0.0; 

T[ 4  *  2 ]  :=  b4  *  El  *  c2  /  sqr(lth); 

T[ 4  *  3  J  :=  b4  *  El  *  c3  /  1th; 

T[ 4  * 4  ]  :=  T[ 3,3  ]  5 
T [ 4  *  5 ]  :=  T [ 2 , 3  ]  ; 

T[ 4  *  6 ]  :=  0.0; 

T[ 5  » 1 ]  :=  0.0; 

T [ 5 » 2 ]  :=  b4  *  El  *  cl  /  (sqr(lth)  *  1th); 
T[5,3j  :=  4  *  2 ]  ; 

T[ 5  » 4 ]  :=  T[ 3,2]  ; 

T[5,5]  :=  T [ 2 , 2  ]  ; 

T[ 5  *  6  ]  :=  0.0; 

T[6,1.  :=  -mu  *  1th  *  sqr(w)  *  ein(th)  /  th ; 

T [ 6 , 2 ]  :=  0.0; 

T[6,3]  :=  0.0; 

T[ 6 , 4 .  • =  0.0; 

T[ 6 , 5 ]  :=  0.0; 

T[6 , 6 ]  :=  cos ( th) ; 

END; 

FOR  i  :=  1  to  m  DO  (*  form  the  matrix  STBGTG  *) 
BEGIN 

FOR  j  :=  1  to  m  DO 
BEGIN 


STBGTG [ 

i 

> 

18 

+ 

j  J 

= 

ID[i,j] 

STBGTG[ 

i 

» 

24 

+ 

j] 

= 

B1  2  [  i  ,  j  ] 

STBGTG [  6 

+ 

i. 

12 

+ 

j] 

= 

T[  i  ,  j  ] 

STBGTG [  6 

+ 

i, 

24 

+ 

j] 

= 

B22 [ i  ,  j  ] 

STBGTG [12 

+ 

i, 

24 

+ 

j  ] 

= 

B32[i,  j] 

STBGTG [ 1 2 

+ 

i, 

30 

+ 

j] 

= 

B33[i,  j] 

STBGTG [ 18 

+ 

1, 

18 

+ 

J  ] 

= 

B35[i,  j] 

STBGTG [ 24 

+ 

i, 

18 

+ 

j] 

= 

T[i,j] 

STBGTG [ 24 

+ 

i, 

24 

+ 

j] 

= 

B46[i  ,  j] 

STBGTG [ 30 

+ 

i. 

24 

+ 

j] 

= 

B32[i,  j] 

STBGTG[30 

+ 

i , 

30 

+ 

j] 

= 

B33[i, j] 

STBGTG [ 36 

+ 

i  » 

18 

+ 

j] 

= 

B35 [ i *  0  ] 

STBGTG [ 84 

+ 

i, 

6 

+ 

0  ] 

= 

T[i,j] 

STBGTG [84 

+ 

i, 

24 

j] 

= 

B1 1  2  [  i  ,  j  ] 

STBGTG [ 90 
FOR  k  :=  1 

+ 

i  , 
to 

18  + 
m  DO 

j] 

= 

T[i,  j] 

BEGIN 

STBGTG [12  +  i,  12  +  j]  :=  STBGTG [ 12  +  i,  12  +  j] 

+  B3l[i,k]  *  T [ k , j ]  ; 

STBGTG [ 1 8  +  i,  6  +  j ]  : =  STBGTG [ 18  +  i,  6  +  j  ] 

+  B34[i*k]  *  T[k, j]  ; 
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STBGTG [ 1 8  + 

i, 

12  + 

S1BGTG [18  + 

i, 

24  + 

STBGTG [ 24  + 

i, 

6  + 

STBGTG [ 90  + 

i, 

12  + 

END  ; 

STBGTG [ 1 6  +  i. 

24 

+  j] 

STBGTG [ 30  +  i. 

6 

+  j] 

STBGTG [ 36  +  i. 

j  ] 

STBGTG [ 36  +  i, 

12 

+  j] 

STBGTG [ 36  +  i. 

24 

+  j] 

ENE 

»; 

END; 

FOR  i 

:=  1  to  m  DO 

BEGIN 

FOR  j  :=  1  to  p  DO 
BEGIN 

STBGTG [ 42  +  i , j ] 
STBGTG [ 48  +  i  ,  j  ] 
STBGTG [54  +  i »  j ] 
STBGTG[60  +  i,j] 
STBGTG [ 66  +  i,j] 
STBGTG [72  +  i,j] 
STBGTG [78  +  i  ,  3  ] 

END; 

END; 


j]  :=  STBGTG [ 1 8  +  i,  12  +  j] 
+  B42[i,k]  *  T[k,j] 

j]  :=  STBGTG [ 18  +  i,  24  +  j] 
+  B4l[i,k]  *  T[k,j] 
j]  :=  STBGTG [ 24  +  i,  6  +  j] 

+  B45 [ i » k j  *  T[ k , j  ] 

j]  :=  STBGTG [ 90  +  i,  12  +  j] 
+  B122[i,kJ  *  T[ k , j ] 

:=  STBGTG [ 1 8  +  i ,  24  +  j] 

+  B36 [ i » j  ]  ; 

:=  STBGTG [ 1 2  +  i  ,  12  +  j  ]  ; 

:=  STBGTG [ 1 6  +  i,  6  +  j]  ; 

:=  STBGTG [ 1 8  +  i,  12  +  j]  ; 

:=  STBGTG [ 18  +  i,  24  +  j  ]  ; 


= 

STBGTG 

[24 

+ 

i»  j] 

= 

STBGTG 

[30 

+ 

i»  j] 

= 

STBGTG 

[36 

+ 

i.  j] 

= 

STBGTG 

[24 

•f 

i,  j] 

= 

STBGTG 

30 

+ 

3  ] 

= 

STBGTG 

:36 

+ 

i*  j] 

= 

STBGTG 

[24 

+ 

i»  j  ] 

FOR  i  :=  1  to  n  DO  (*  initialize  the  matrix  DIAG  *) 
BEGIN 

FOR  j  :=  -ml  to  m2  DO 
BEGIN 

DIAG [ i  ,  j  ]  :=  0.0; 

END; 

END; 


FOR  i  :=  0  to  95  DO  (*  form  the  matrix  DIAG  *) 

BEGIN 

FOR  j  :=  1  to  36  DO 
BEGIN 

DIAG [ i  +  1,  -19  -  (i  mod  6)  +  j]  :=  STBGTG [ i  +  1,  ; 
END; 

END; 


wri te  ( Ls t , w) ; 

bandet ( n  ,m1  , m2 , DI AG , det) ;  (* 

writeln  (Lst,det);  (* 

:  *  w  +  1 0 . 0  5  (  * 


call  procedure  bandet, 
write  det  , 

and  add  an  increment  to  w 


w 


